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The coherent potential approximation has historically allowed the efficient study of disorder effects 
over a variety of solid state systems. Its original formulation is however limited to a single-site or 
uncorrelated model of local substitutions. This neglects the effects of correlation and short range 
ordering, often found in realistic materials. Recent theoretical work has shown how to systematically 
address such shortcomings, for simple materials with only one element per unit cell. We briefly 
review the basic ideas of these developments within the framework of multiple scattering theory, 
and suggest their generalization to materials with complex lattices and possibly different types of 
disorder. We illustrate this extension with an example of local environment effects in the exotic 
Hapkeite Fe — Fe^ .5%5«62.5% compound. 

PACS numbers: 71.23.-k, 71. 15. Mb, 71.15.Ap, 71.15.Dx, 71.20.Be, 71.20.Gj, 71.20.Lp, 75.50. Bb, 75.30.Hx, 
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I. INTRODUCTION 

The first-principles treatment of solid state systems in 
the presence of disorder poses an extra degree of tech- 
nical and numerical difficulties, due to loss of crystalline 
periodicity.^ Some deviations from such ideal scenario are 
however to some extent always present: real samples al- 
ways contain various kinds of defects, which may have a 
significant influence on the fine details of their electronic 
structure. 

Moreover, the general framework of a substitutional 
model of disorder, where occupation of any lattice site 
by alternative species is only determined in probabilis- 
tic terms, applies to a variety of other physical problems 
of great practical interest. In addition to the study of 
metallic alloys J*"^, this includes among various examples 
doped semiconductors, but also random spin systems^ 
such as magnetic materials above a critical ordering tem- 
perature^ and further scenarios. 

One efficient technique to tackle these problems is pro- 
vided by the coherent potential approximation (CPA) 8 
in its multiple scatteri ng Korring a, Kohn and Rostoker 
(KKR) implementation B l 1 IE 1 1 1 1 12 1 The theory operates by 
examining the actual physical system in terms of a com- 
putationally more amendable effective medium, which is 
determined self-consistently for each energy so as to best 
incorporate average disorder effects. It has been shown 
to satisfy a strict set of theoretical requirements, re lated 
to convergence and analyticity of a proper solution) 13 ! 14 ! 
to recover exact results in all the appropriate limits^ 
and to provide qualitative, systematic improvements over 
cruder models of disorder-^ at a reasonable increase in 
computational demands. 

On the other hand, reliance on concentration alone to 
model the average presence of different chemical or mag- 
netic species cannot account for possible higher order as- 
pects of their actual distribution, as found in a real mate- 
rial!^ Local atomic substitutions in fact typically extend 



their effect up to an extended length scale, which may or 
may not correspond to that of crystalline periodicity, or 
other forms of long-range order (LRO) in a sample. 

Attempts to improve the original single-site CPA for- 
mulation to also describe such short-range order (SRO) 
or intermediate disorder scenarios have prompted a large 
body of theoretical efforts.^ On the one side can be 
placed developments to include the case of complex unit 
cell materials, where some forms of disorder coexist with 
undisturbed periodicity on the other sublatticesP^ From 
a complementary point of view, there have been also at- 
tempts to explicitly re-insert local higher order correc- 
tions in the original theory for simple systems, on top of 
the previous average -only results.^ 

Along the latter line of research, the non-local co- 
herent potential approximation (NLCPA) has been re- 
cently developed based on insights provided by the Dy- 
namical Cluster Approximation (DCA)P3 The method 
has shown its capability to des cribe SRO effects both 
within model Hamiltonian studies^2iH2ll ailc j when reim- 
plemented within a completely self-consistent DFT-KKR 
framework for realistic systemspSHUl 

This technique inherits however from its derivation a 
practical restriction to simple, one atom per unit cell lat- 
tices in Strukturberichte Ah, A2 or Al structures. In this 
paper we reexamine the foundations of such approach, 
and proceed to generalise it to also include the case of 
complex unit cell, multiple sublattices scenarios, in ar- 
bitrary geometries. Following previous discussion for a 
model Hamiltonian application} 2 ^ this work considers in 
particular the problem from the prospective of multiple 
scattering calculations for actual materials. Many poten- 
tial examples spring to mind including high T c supercon- 
ductors, multiferroics, candidate compounds for hydro- 
gen storage etc. To illustrate some of the effects which 
the new method can describe, we presently apply it how- 
ever to the specific case of the Hapkeite iron-silicon min- 
era j29jt3iJ Thig ahoy i s naturally formed in lunar soil by 
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impact-induced, vapor phase deposition, and turns out 
to be a particularly illustrative example for current dis- 
cussion, thanks to its particular lattice simplicity. This 
follows a Pm3m symmetry, where one of the two dis- 
tinct crystallographic positions is periodically occupied 
by Fe atoms, while the other can randomly host cither 
iron or silicon. This kind of intermixing highlights the 
sensitivity of the transition metal not only to the overall 
concentration of the Si atoms, but also to specific local 
environment effects. 

The paper is organized as follows. The essential as- 
pects of the original CPA method are briefly reviewed, 
and its single-site limitations are tracked down in par- 
ticular to the step where the conditional averages over a 
variety of local realizations of lattice occupations are eval- 
uated. Our analysis points out how higher order correc- 
tions, already quantified by Mills et al. for model Hamil- 
tonians in the travelling cluster approximation (TCA) P=l 
are missed out in the original theory. We recap from a 
similar angle the alternative strategy to reinclude these 
adopted by the NLCPA (section pi)); and the independent 
generalisation of the single-site CPA, proposed instead 
by Pindor et al. 19 to describe complex unit cell cases 



(MS-CPA, section II C I. Finally, we present our unified 



solution for the treatment of LRO and SRO scenarios 
in complex lattice materials as a formal merger of the 
MS-CPA and NLCPA (section [In} . 

This multi-sublattice, non-local generalization of the 
theory (MS-NLCPA) is first put to the test by artificially 
examining the A2 (Al) disordered Cu\- c Zn c metallic al- 
loy as an example of a multi-sublattice B2 (LI2) material, 
and showcasing the full equivalence of our MS- NLCPA 
results to those of the original 

NLCPAEHl (section |hTa). 



We then demonstrate the new functionality of the 
method with an application to the Hapkeite Fe — 
Fe c Si\- c compound. We explore here in particular the 
connection between electronic and magnetic properties 
of the iron atoms as a function of their average local 
environment, under a fixed experimental concentration 
c = 37.5%. 

The main aspects of these developments are finally 
summarized in section HVl 



II. 



EFFECTIVE MEDIUM STUDIES OF 
DISORDERED SYSTEMS 



We briefly review the basic ideas of the theory by re- 
ferring to a generic binary alloy in the form: A\- C B C . 
Our interest lies with observables that involve a statis- 
tically significant number of randomly occupied lattice 
sites, over which averaged properties are assessed. Typi- 
cal experimental techniques include those of photoemis- 
sion and nuclear magnetic resonance spectroscopies; elec- 
trical and optical transport; estimates of magnetic and 
compositional transition temperatures; and so on. 

In all these cases, it is a great many configurations of 
atomic species A and B within different local arrange- 



ments that give rise to the aggregated, observable result. 

We approach the computational challenge posed by 
these scenarios by replacing the original problem with 
an effective medium description. This represents a fur- 
ther step of abstraction in the context already established 
by density functional theory, in which one-electron -like 
Kohn-Sham orbitals are evaluated starting from an ap- 
proximation of the actual many-body potential gener- 
ated by all charges, in the specific arrangement of various 
atomic sites on a lattice. 

In a multiple scattering solution, these "scatterers" are 
decoupled in the effect produced by local interactions, 
and the global spatial distribution. Such separation is 
used to efficiently determine the electronic propagator 
Q(k,e), and from it the electronic charge density. This 
provides the necessary input for a further iteration of the 
Kohn-Sham scheme, until convergence is achieved and 
additional observables may be determined^. 

When some form of disorder disrupts periodicity of oc- 
cupation over the various lattice sites, the CPA applies 
to the above general scheme an additional demand that 
the system Green function should reproduce the average 
properties of the actual material, evaluated over a large 
enough portion of the bulk. We briefly review the basics 
of this technique to also highlight one of its specific lim- 
itations, before examining some suggested amendments. 



A. The coherent potential approximation and its 
single-site nature 

In the original derivation by Faulkner and Stocks,^ 
an approximation for the ensemble averaged propagator 
G(k, e) is offered in terms of an effective medium scatter- 
ing path operator r(e). This mathematical object con- 
tains all specific properties for the system, and is effi- 
ciently and conveniently represented in a spherical har- 
monics and lattice site basis, denoted here in its matrix 
form over the L = (£, m) and site i indices as r(e). 

The CPA provides a scheme to obtain such quan- 
tity from the condition of no extra scattering coming 
on the average from any embedded impurities. In prac- 
tice this constraint can be enforced starting for each en- 
ergy e from the T-matrix approximation (ATA) Ansat J^, 
which posits for the t-matrix associated with every site 
the average form: 



a 



AovB^i(e) = (l-c)t A (e) + ct B (e) (1) 



If we assume that the CPA effective medium should be 
on average the same on each lattice position n (as often 
indicated in the literature through the convention n = 0, 
here and in the following), we can adopt a Brillouin zone 
integral representation and restrict ourselves to the finite 
site-diagonal part of r (e) : 



J dk(m(e)-G{k,e] 



(2) 
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where fi is the Brillouin zone volume, m(e) = \ (e), and 
G(fc, e) represents the structure constants matrixPS 

A self-consistent prescription can then be set up, by 
also requiring that a concentration weighted sum of the 
conditional averages over constituents in direct space 
should also lead to the same result, 



r n '"(e) = (1 



c)r n /\e) + cr^ n (e) = Y / C a r n a n (e) (3) 



where T^ ,n (e) describes an impurity of type a, embed- 
ded in the CPA medium at site n. This quantity can be 
obtained through application of the corresponding pro- 
jectorPl 

D a (e) = (l + r"' n (e)(m Q ( e )-m( £ ))) _1 (4) 

for m Q (e) = ^(e), so that: 

^ n (e) = D a (e)r^(e) (5) 

and iteration between eq. [2] and [3] until both pre- 
scriptions converge to the same scattering path opera- 
tor T™'™(e) provides the desired CPA description of any 
average site, within the supposedly infinite bulk. 

We now follow Mills et aiPSI in calling attention to a 
basic assumption of this approach. The conditional en- 
semble average of N — 1 randomly occupied sites has been 
so far worked out by simple factorization into N fully un- 
corrected substitutions. However when more than one 
disordered contributions are to be considered together, a 
proper averaging procedure should also contain higher or- 
der terms. In the generic example of N random variables 
X\. X 2 , . . . , Xn, these corrections can be rec ursively de- 
fined in the form of cumulant averages (CA)|i21M3Z| 



<*i) 
(X 1 X 2 ) 

(XiX 2 X 3 ) 



= (X 1 )(X 2 ) + (X 1 X 2 ) C 
= (X 1 )(X 2 )(X 3 

+ (X 1 X 2 ) C (X 3 ) C + (X 2 X. 

+ (X 1 X 3 ) C (X 2 ) C 



(XiX 2 X 3 ) - 



{XtX 2 



.X 



N 



(X 1 )(X 2 )(...)(X N ) + 



(6) 

which are rigorously null only in a single site case, and 
otherwise given by (XiX 2 ) c = {X X X 2 ) - (X l ){X 2 ), and 
so forth. 

In the present context, these terms reflect the influence 
that placement of an impurity a on site n exercises on its 
neighborsP^In a SRO regime, such an impact decays over 
a smaller distance than the size of the bulk, and does not 
correspond to the effective medium periodicity, implicitly 
assumed in eq. [2] Significant effects can however arise, 
which are missed in the cruder single-site factorization of 
the original CPA theory. 

Our approach to overcome such limitation picks up 
from the non-local extension of the methodpS which we 
now briefly proceed to review. 



B. The Non-Local CPA solution 

The NLCPA operates as well within the basis of a fac- 
torized evaluation of the effective medium, combining the 
contributions from different substitutional species a. 

In this formalism, however, these are now considered 
over clusters of N c > 1 sites, occupied by sets 7 = 
a 2 , ■ ■ • , CKjv c } °f possibly correlated impurities. A 
multi-site probability distribution P{"f) generalizes then 
the role of the single-site concentration c Q , and can ac- 
count for different forms of SRO over a length scale fixed 
by the size of the cluster. When N c = 1, the method 
reverts to the CPA treatment, from which it also inherits 
all the analytical features of a Herglotz solution. It be- 
comes however exact in the limit of cluster size extending 
to the whole bulk, N c = N — > 00. 




FIG. 1: Schematic depiction of a multi-site cavity with N c = 4 
elements (darker gray). The cavity is evaluated at an arbi- 
trary position from the effective medium (lighter gray). Vari- 
ous substitutional impurities are placed in different, possibly 
correlated arrangements 7, according to a multi-site proba- 
bility distribution P(j)- The different labelling of distinct 
site occupation pertains to a NLCPA (cti index, N su b = 1, 
iV c = 4) or MS-NLCPA (ai, s index, N sub = 4, N c = 1) de- 
scriptions, which can be equivalently deployed and lead to the 
same results in ad hoc degenerate test cases (see text). 

Technically, this is accomplished through a modified 
version of eq. [2] and [3j based on an enlarged but finite 
matrix structure for the scattering path operator r(e). 
Here a second underline denotes now a blockwise exten- 
sion of the original CPA matrix into an extra pair of in- 
dices, labeling different I,J— 1,...,N C cluster elements 
in the bulk (Fig.[T]). The corresponding impurity projec- 
tors become: 



(l + f(e)( 



m : 7 ( e ) 



(7) 



4 



for m (e) = t 1 (e) a matrix diagonal in cluster site in- 
dices, describing the single site scattering from atoms at 
sites I within the cluster, according to each particular 
configuration 7: 



r y (e)=£ 7 (e)£(e) 
Similarly to eq. [3J we then assume: 

£(e)=^P( 7 )r (e) 



(8) 



(9) 



or may equivalently express the contribution of each con- 
figuration 7 in terms of a cavity scattering path operator 
r ca "(e), such that P 



Z cav (e) = ||(e) -£ _1 (e) 
r 7 (e) = (^(e)- S c -(e))' 



(10) 



Care is needed when extending eq. [2] in reciprocal 
space, to correctly describe a finite cluster still repre- 
sentative for the whole bulk. Jarrell et ai.p^ and sub- 
sequently Rowlands et ai.p' have shown how this can 
be accomplished through a partitioning of the original 
Brillouin zone domain of integration respectful of the un- 
derlying lattice symmetries. This is now coarse-grained 
into N c tiles around a discrete set of cluster momenta 
K n , which remain defined onl y u p to an arbitrary phase 
factor, further discussed below"*. 

The corresponding scattering path operator is hence 
obtained in the modified form: 



r{K n ,e) = ^J dk(t \e)-SG(K n ,e)-G(k,e) 



L / dk(m(K n ,e)-G(k,e) 



(11) 

where stands for the integration volume of each of 
the N c tiles, defined so as to preserve on-average transla- 
tional invariancepS SRO corrections 8G_(K n ,e) for each 
cluster momenta are made to appear, and can be col- 
lected into the now K n -dependent m(K n , e) construct 

of cq. [IT] together with the original \ 1 (e) . 

In practice, the algorithm begins again from a block- 
wise expression of the original ATA Ansatz of eq. [T] 
initially without these extra terms. It iterates between 
the results of eq. [9] and eq. fTT] through the coupled set of 
lattice Fourier transforms!^ 



(12) 



until satisfactory convergence to the same effective 
medium description is achieved. Here Rj denotes the 
position of a site within the embedded cluster, and the 
Brillouin zone tile centers K n are related to these vectors 
such that i J2 n e lK "- {R '- R ^ = 5 LJ . 

The tiling procedure places however some constraints 
on the allowed sets of such N c calculation parameters.^ 



To this date the technique has been developed only for 
simple lattices with just one element per unit cell, in the 
case of Strukturberichte A^, A2, Al geometries Fur- 
thermore, this coarse-graining of the Brillouin zone leads 
to discontinuities in the k -dependent integ rand of Eqfl2| 
whenever tile boundaries are crossed 1^1221 \\fe note here 
however a recent suggestion to overcome such limitation 
through an additional averaging step over various phase 
choices, to reobtain a smooth k -dependence^. 

Before proceeding to develop the multi-atom per unit 
cell generalisation of this KKR-NLCPA approach, we 
shall now briefly recall for clarity also the complementary 
starting point of the single-site CPA applied to complex 
lattices. 



C. The multi-sublattice CPA for complex lattices 
in arbitrary geometries 

We consider here the extension of the single-site CPA 
proposed by Pindor et al. for multiple sublattices com- 
pounds (MS-CPA) .121 This theory is designed to cover 
cases of "periodic disorder" , where a material with s = 
1,2,..., N su b crystallographic positions in the unit cell 
may randomly host more than one atomic species per 
site. A prototype N su f, = 2 example can be given by the 
ax = A or B, cti = C or D abstract compound, in the 
general formula: (A 1 _ Cl ,B Cl ) - (Ci_ C2 , D C2 ). 

In this case, the reciprocal space prescription for the 
site -diagonal scattering path operator of eq.J^becomesP^l 

r SiS ,(e) = ^/dfc(r 1 ( e )-G(fc, e ))~ 1 *„/ = 

Y~ / s,s' 

= j l fdk(m(e)-G(k 7 e)y 1 S ss/ =f s (e) 

\ / s,s' 

(13) 

where t (e) is the block-diagonal matrix with elements 
t s S SjS i in sub-lattice space, and now G(fc, e) describes the 
structure constants for the free electron propagation in 
the complex lattice. The self-consistent determination 
of the effective medium starts once again from the ATA 
Ansatz of eq. [Tl modified to have: ^(e) = (1 — ci)t A {e) + 
Cit B (e), andtj(e) = (1— c 2 )t c (e)-|-C2i D (e), in this generic 
example of a binary alloy with two sublattices. 

This prescription is complemented by considering in 
direct space 



t 1 (e) = (l-c 1 )r 1 A (e) + c 1 r 1 B (e) 
r 2 (e) = (1- c 2 )t 2 c (e) + c 2 T 2 D (e) 



(14) 
(15) 



with contributions from the different impurities that may 
be found on each sublattice now given by the projectors: 

D as (e)=(l+r(e)(m as (e)-m s (e))y 1 (16) 
so that in general: 



r s a (e)=D a (e)r s (e) 



(17) 
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We note that the formulation remains sub-lattice 
block-diagonal at all steps, despite the complete N su b x 
Ngui, block matrix inversion appearing in eq. [13] When 
compared with the NLCPA of cq. |9"fTT] it can be seen 
to describe at most LRO or periodic disorder cases, but 
no local effects from SRO including those linking the oc- 
cupancies of the different sub-lattices. As expected, this 
development also reverts back naturally to the single-site 
CPA for N su b = 1. It shares with the original formula- 
tion also lack of particular restrictions to specific lattice 
geometries, another element that we wish to retain in our 
generalisation proposal. 



III. MERGING THE TWO TREATMENTS 

Our strategy for combining respective benefits from the 
two theories of sections |II B| and |II C| will be again based 
on a self-consistent procedure to track multiple length 
scales in an extended cavity, as depicted for the NLCPA 
in figure [I] We still resort in particular to considering 
multi-site substitutions governed by a richer probability 
distribution P(7), suitable to describe different forms of 
short and long range ordering. This general scheme is 
derived from the NLCPA, and indeed the following anal- 
ysis will be based on the formal demand that the new 
generalized theory should recover the results from such 
starting point, in the appropriate limit. 

We intend to also match this extended model of disor- 
der with assumptions similar to those underneath eq. |13| 
where now all terms are carefully preserved and no off- 
diagonal contributions in sublattice space get discarded. 
The resulting richer description of the effective medium 
contains additional non-diagonal t-matrix terms t s s , (e) , 
similar to the |/j(e) contributions in the original NL- 
CPAPUWhen both developments are combined, the scat- 
tering path operator is also extended, to ultimately ob- 
tain a matrix r(e) now labeled in general by angular mo- 
mentum, tile and sublattice indices. 

This allows the formalism to track the effects of possi- 
bly correlated substitutions over N cuto f / = N c x N su f, ele- 
ments' length scales. Special care must however be given 
to the important role now played by the off-diagonal 
blocks of the integrand from eq. [13] 

We consider here for simplicity a N c — 1 example. 
The general expression for the extended scattering path 
operator T ss i (e) can be written as: 



\ / ss' 



(18) 



for m(e) = I (e) at first given by the usual ATA pre- 
scription. 

The term G(fc, e) denotes now however modified struc- 
ture constants matrix blocks G s s ,(k, e), obtained from 
the original one by considering the extra sublattice - 
dependent phase modifier: 



G sa ,(k,e) = G sa ,(k,e)e- ik < r *- r ^ 



which can be postulated by following Banachiewicz's the- 
orem for the blockwise inversion of a square matrix. Here 
r s is the position vector of site s within a unit cell, and 
an intuitive interpretation for the action of such correc- 
tion can be simply illustrated considering the two cases 
of Strukturbericht A2 and B2 lattices. 

In this latter geometry, the off-diagonal structure con- 
stants matrices G^^, l {k, e) are related to the correspond- 
ing ones for a BCC case according tcff^l 

G?f l (k,e) = (G BCC (k,e)-G?f l (k,e)) e +*^-r 2 ) 

(20) 

Applying eq. |19| leads hence to obtain, for the non trivial 
off-diagonal blocks: 



= G^ cc (k,e)-G?f l (k,e) 



(21) 



where r 1 = (0,0,0), r 2 = (1/2,1/2,1/2) are possible 
basis vectors for the complex unit cell of an A2 lattice. 

Such cancellation of the exponential from eq. [20] oper- 
ates in other words by removing the sublattice -diagonal 



free propagation modes G u s l (k,e) 



G 



CsCl 
22 



(k. 



from 



the envelope of a generic expression for G ss ,(k, e), ini- 
tially designed to account for both diagonal and off- 
diagonal hopping processes (Fig. [2]). 



(a) 



N su h — 2 



7 V — ^ 



// (ft) 



G a=al (k,e) = G 3=s ,(k,e) 



G S * S \k,e)^G^\k,e) 



(19) 



FIG. 2: Schematic separation into distinct propagation modes 
for purely sublattice diagonal (s = s', (6)) and off-diagonal 
(s / s', (c)) free electron hopping processes, in a N S ub = 2 
sublattices case (blue and red) . Different length scales can be 
tracked consistently with the comparable NLCPA treatment 
of a N c = 2 simple lattice scenario, by decomposing the ag- 
gregated results of panel (a) into the separate contributions 
of panels (b) and (c). 

For the s = s' terms, there is no exponential argu- 
ment and the original expression of panel (a) remains un- 
changed. The suggested modification of eq. [19] however 
ensures that, when integrating also for the new s =/= s' 
contributions, such blocks of the scattering path opera- 
tor keep resolving well distinct length scales effects which 
may be sensible to sublattice-sublattice couplings. Short 
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range ordering can then be set-up in direct space by 
straightforward combination of eq. 7]|9 and T4][T6 as de- 
tailed below. 

We find thus in general: 

r St AK n ,e) = § J QKn dk (t = (K n , e) - k) 

U, s ,j,Ae) =^xi^sA K n^ +lK -- [R '- R ^ 

(22) 

as a result of applying, if required, the NLCPA lattice 
Fourier transforms of eq. [12] with Rj , Rj and K n now 
referring to superlattice vectors. 

This larger scattering path operator r(e) is again 

determined self-consistently by solving for contribu- 
tions from different extended cavity occupations 7 = 
{ai.i, . . . , ai yS , . . . ,ciN c ,N sub }, through a generalization of 



eq. |7fl6l or the equivalent extension of eq. [TO 

— (e) = (rH^-^ie) 
so that as in eq. [9} 



— *■ — * 



(23) 



(24) 



Eqs. 22]|9 are then iterated self-consistently, until con- 
vergence to the same effective medium description is 
reached, and other observables of interest may be com- 
puted. 



A. Validation tests 

We propose that a practical approach to validating the 
above procedure can be conceived by purposedly exam- 
ining simple one atom per unit cell lattices, taken as par- 
ticular instances of multi-atom per unit cell materials on 
a superlattice. We then require that our MS-NLCPA 
should reproduce the original NLCPA results in this de- 
generate limit, over all equivalent SRO scenarios. To this 
en d we consider for instance a typical Cu 50 %Zn 50 o/ al- 
lojIMH] i n either a BCC (A2) or FCC (Al) phase, and 
explore different local environment regimes with the two 
techniques. 

The set of fully comparable structures are reported in 
table [U A "degenerate" B2 setup, equivalently hosting 
either Cu or Zn atoms on both sublattices (N sub = 2, 
N c = 1), can also be examined in a NLCPA study of a 
A2 unit cell (N su b — 1, N c = 2). Similarly, a LI2 system 
{N sub — 4, N c — 1) can also be set up to mimic the 
equivalent (N sub = 1, N c = 4) Al NLCPA setup. In 
each case the results should be identical. 

This equivalence is confirmed for all the relevant ef- 
fective medium quantities throughout the self-consistent 
calculation. Following Rowlands et aZp3, we show here 
in particular the density of states for three very different 



Computational budget: 
(Nr. of lengthscales) 


Comparable descriptions: 
Lattice type: N c N su b 


2 


BCC or A2 2 1 
CsCl or 52 1 2 


4 


FCC or Al 4 1 
CuzAu or LI2 1 4 



TABLE I: Comparable structural decompositions of the same 
two physical systems (first and second row). When each sub- 
lattice is constrained to host the same elements with the same 
probabilistic distribution and over matching lattice parame- 
ters, the NLCPA and MS-NLCPA are expected to coincide at 
the level of all effective medium quantities. 



SRO prescriptions for the two physical cases described in 
columns 1 and 2 of table [I] (Fig. [| and [5} . 

In all cases, the SRO = regime of a fully 
uncorrelated probability distribution pSRO=0(^ _ 
Il/=i ris^i 6 c ( a /,s(7)) corresponds to a factorized eval- 
uation of the various elements' concentrations c through- 
out the cavity, for a/^f/y) the atomic species found on 
site /, s. As already observed in the NLCPA develop- 
ment, this model of disorder shows in general little dif- 
ference from the outcome of straightforward single-site 
CPA calculations. 

Additionally however, the extremal cases of full, local 
clustering (SRO = +1): 



P 



SB.O=+l 



(7) = c(a(7)) if oti lS — a V/, s 



(25) 



or full local ordering (SRO = — 1; details in figures 
captions) may also be successfully compared, as well as 
any other types of local environmentspH The new MS- 
NLCPA appears to incorporate the NLCPA results as 
desired. 



B. Application to the B2 Fe — Fe 37 5 %Sig2.5% 
Hapkeite mineral 

As more realistic application of the new method, we 
examine now the example of the Hapkeite phase of the 
Fe - Si systenPHSl. 

Our choice of such compound is partly based on the de- 
sire for a relatively simple yet representative case study 
to demonstrate the new capabilities of the formalism. 
The above mineral is characterized by a space group 
Pm3m (aiat — 2.81A) closely packed structure^, with 
a N su b = 2 complex unit cell hosting Fe on the first sub- 
lattice ri = (0,0,0), and cither Fe or Si in the body 
central position r 2 = (0.5,0.5,0.5). 

We focus here in particular on the effect of SRO, for 
a fixed experimental concentration ratio of Fe : Si — 
3 : 5 on the disordered sublattice, the other sublattice 
being regularly occupied by iron atoms. Three scenarios 
of completely uncorrelated disorder (SRO — 0), short- 
range clustering (SRO = +1) and local ordering (SRO = 
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FIG. 3: The density of states versus energy (color online) 
for a Cu 50 %Zn 50 % alloy in both a BCC and CsCl struc- 
tures, evaluated for N su b x N c = 2 in the 3 extremal short- 
range ordering regimes of no correlation (SRO = 0, in blue: 
P(7i = {Cm, Cm}) = 25%, P( 72 = {Zti.Cm}) = 25%, 
P(f3 = {Cu,Zn}) = 25% and P(7 4 = {Zn,Zn}) = 
25%); complete bias towards alike neighbours (clustering, 
SRO = +1, in red: P( 7 i = {Cm, Cm}) = 50% and P(~f 4 = 
{Zn,Zn}) = 50%); complete bias towards unlike neighbours 
(ordering, SRO = -I, in green: P(~f 2 = {Zn,Cu}) = 50% 
and P(7 3 = {Cu,Zn}) = 50%). The NLCPA results are the 
dotted lines whereas the continuous line shows those of the 
MS-NLCPA. In all 3 cases, the MS-NLCPA CsCl results are 
indistinguishable from the BCC NLCPA ones. 



FIG. 4: The density of states (color online) for a Cu5o%Zn 5 o% 
alloy in both a FCC and CM3AM (LI2) structure, evalu- 
ated for N 3U b x N c — 4 in the 3 extremal short-range or- 
dering regimes of no correlation (SRO = 0, in blue: P(~fi = 
{Cm, Cm, Cm, Cm}) = P( 72 = {Cm, Cm, Cm, Zn}) = ... = 
P(7i6 = {Zn, Zn, Zn, Zn}) — 6.25%); complete bias towards 
alike neighbours (clustering, SRO = +1, in red: P(fi = 
{Cm, Cm, Cm, Cm}) = 50% and P(~f 2 = {Zn, Zn, Zn, Zn}) = 
50%); complete bias towards unlike neighbours (ordering, 
SRO = -1, in green: P( 7 i = {Cu, Zn,Cu, Zn}) = 50, 
P(72 = {Zn,Cu,Zn,Cu}) = 50%). The NLCPA results are 
the dotted lines whereas the continuous line shows those of 
the MS-NLCPA. In all 3 cases, the MS-NLCPA Cm 3 .4m results 
are indistinguishable from the FCC NLCPA ones. 



— 1) are again explored, similarly to the validation studies 
discussed in Fig. [3] and |4j 

We set the maximal correlation length on which 
such possibilities are considered to the size of a multi- 
site cavity Nf^ s to ^ = 8. This setup can be eas- 
ily obtained in direct space as the 2x2x2 super- 
cell of the original B2 structure. In a fully clus- 
tering scenario, the N to t = 2 8 = 256 local config- 
urations 71 = {Fe s=1 , . . . , Fe s=8 , a s=9 , . . . , a s=1% } for 
a s = Fe or Si, that specify all the possible fill- 
ings, are restricted to just N^f 0=+1 = 2 cases: 
7 SRO=+i = {Fe s=1 ,. . . , Fe s=s ,Fe s = 9 , Fe s = 16 } with 
probability P( 7l si? ° =+1 ) = 0.375; and 7| fl ° =+1 = 
{Fe^ 1 ,...,^^ 8 ,^^ 9 ,...,^^ 16 } with probability 
P(l2 R ° =+1 ) = °- 625 - Alternatively, full local order- 
ing can be defined as the restriction to those configu- 
rations, where strictly as many Si atoms as mandated 
by the concentration can occupy the central positions, 
and all the other ones are hosting iron. The desired ele- 
ments' ratio is then recovered by giving to each of such 

Ntot 0-1 = ( |j ^ =56 arrangements a uniform proba- 
bility distribution: P( 7 f fl ° =_1 ) = . . . = P(7 5 S 6 RO=_1 ) = 
l/Nf ot °=~ l ~ 0.018. Finally, the case of no restrictions 
corresponds to the fully uncorrelated, SRO — regime, 
which is implicitly assumed in the original MS-CPA. 

In all these scenarios, a full deployment of the formal- 
isnPH in the KKR-CPA scheme generates a new charge 
distribution associated with each substitutional configu- 
ration 7j, through iteration until convergence of a stan- 



dard Kohn-Sham LDA-DFT scheme. This allows to com- 
pute element -resolved magnetic moment variations as- 
sociated with disorder and SRO, as well as other observ- 
ables such as in this case the spin-resolved DOS. These 
results are shown in Fig. [5| 

We observe most significant variations for the cluster- 
ing scenario (SRO — +1, Fig.[7|), as opposed to the other 
two. The most severe constraint on local microscopic 
compound realizations (from Nt t°~° — 256 down to 
N tot'° =+1 = 2 alternatives) translates into DOS plots 
with a more defined outlook, in terms of sharper fea- 
tures across the whole energy spectrum. This can be 
easily understood noting how the limiting case of com- 
plete periodicity or ideal LRO would in fact appear when 
deploying the method for a N to t — 1 setup, thereby fully 
suppressing disorder. Conversely, the local ordering case 
SRO = — 1 appears to be quite sufficient in capturing 
the essence of electronic structure contributions from the 
two differently occupied sublattices. Results from an un- 
correlated model of disorder (SRO = 0, Fig. [6]) are in 
fact closely mimicked by evaluating the 56 contributions 
of variously permuted Fe and Si elements in the appro- 
priate lattice sites for this locally ordered choice. 

The study can lead however to further insight as well 
into the local environment-resolved impact of distinct 
constituents. Considering for example the atomic mag- 
netic moment, the fully uncorrelated MS-CPA case pro- 
vides reference values of 0.48 [hb] for the iron atoms 
F e ord~ >= ° on t ne completely -Fe-filled sublattice, whereas 
those hosted on the other sublattice in competition with 
silicon atoms, Fe^ o=0 , have a larger moment of 2.65 



1Mb]- 

These properties change quite markedly when clus- 
tering is instead imposed (SRO = +1, Fig. [7|. We 
note that, while the magnetic moment of the disor- 
dered sites remains substantially unchanged (2.64 [/is] 
for Fe^° =+1 ), the Fe^ d ° =+1 atoms on the other sub- 
lattice acquire a magnetic moment of 1.09 [(J-b] when 
surrounded again by Fe (configuration 7f fl °~ +1 ), but 
only 0.19 [/j-b] hi a .Si-only local environment (configu- 
ration 72 )■ A strong hybridization between elec- 
tronic levels in the former case leads in fact to the en- 
hancement of local magnetic moment of the periodically 
repeated Fe^ d °~ +1 atoms, as highlighted in the corre- 
sponding spin -resolved DOS shown in Fig. [5] 



Spin -rcsovcd DOS N clltoft =8 Fe-Fe 37 ,5%,Sig2.5% 




-9 -8 -7 -6 -5 -4 -3-2-10 1 2 



E - Ep [eV] 

FIG. 5: Density of states for the Hapkeite mineral Fe — 
Fe37.55i62.5 in the B2 phase, with one sublattice periodically 
occupied by Fe, the second one hosting either Fe or Si un- 
der different SRO assumptions (color online). We contrast 
in particular a fully uncorrelated (SRO—O, red) probability 
distribution of constituents with the case of local ordering 
(SRO=-l, green) or local clustering (SRO=+l, blue). 



Further differences arc evident in the SRO = — 1 sce- 
nario. A similar self-consistent calculation attributes to 
Fe^° __1 atoms a magnetic moment of 0.58, 0.52 or 
0.48 [/is], emerging as clearly distinct values across the 
56 local configurations that are now taken into account. 
This can be correlated with the corresponding realization 
of three specific types of local geometries. The highest 
magnetic moment is for instance associated with those 
8 cases where complete iron layers along the [111] direc- 
tion (and equivalent ones such as [1,1,1] etc.) are formed 
within the 2x2x2 supercell. The lowest magnetic mo- 
ment is found instead for local arrangements realizing 
any of the 6 x 4 = 24 (1, 1,0) (and permutations) iron 
layers, which also include a Sip e defect; and a similar 
analysis can also be repeated for the intermediate value 
of 0.52 [hb\- For the non-periodic Fe d ^ 0= ~ 1 atoms, the 
magnetic moment appears also close to the SRO = 
value, with less significant local environment-dependent 
variations around an average value of 2.69 \jjlb\- 



Element and spin -resolved DOS Fe-Fc 37 ,5%,Si 62 .5%, SRO=0 
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FIG. 6: (color online) Atomic resolved densities of states and 
local magnetic moments for the Fe — Fe 37 . 5 o /o Si 6 2.s% alloy in 
the fully uncorrelated case (SRO = 0): sublattice periodically 
occupied by Fe (a); sublattice hosting either by Fe (b) or Si 
(c) atoms. 



IV. CONCLUSIONS 

Most materials are affected by some degree of disorder 
and in many crystalline systems this affects differently 
a subset of the sublattices. Often properties nominally 
associated with the chemical composition of one sublat- 
tice can be modified by the extent and nature of disorder 
on another. This effect and more general cases of short 
range ordering regimes can be overlooked by the origi- 
nal CPA, due to the single-site nature of the theory. In 
this paper we briefly reexamine further extension to this 
general method, and suggest a combined approach to ef- 
ficiently perform first-principles studies on this class of 
systems. 

Previous results for a reference Cu^%, Zn^% alloy in 
a A2 (Al) I B2 (Ll 2 ) structure could be used to pro- 
vide strict validation of the new formalism. Further- 
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more, this could be put to test in its new capabilities, 
considering here the particularly simple example of a 
Fe — Fe 37 5 %Si e2 .5% Hapkeite compound examined from 
the point of view of basic electronic and magnetic effects 
induced by different forms of SRO. 
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FIG. 7: (color online) Atomic resolved densities of states and 
local magnetic moments for the Fe — Fez-j.sSim.b alloy in 
the case of full clustering (SRO = +1, blue): sublattice pe- 
riodically occupied by Fe in j^ro^+i ^ and ^sro=+i ^ 
configurations; sublattice hosting either by Fe (c) or Si (d) 
atoms. 



